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Abstract. We prove that the inverse spectral mapping reconstructing the square 
. integrable potentials on [0, 1] of Dirac operators in the AKNS form from their spec- 

CNJ ' tral data (two spectra or one spectrum and the corresponding norming constants) is 

, analytic and uniformly stable in a certain sense. 
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1. Introduction 

The main goal of this paper is to estabhsh analyticity and uniform continuity of 
solutions in the inverse spectral problems for a certain class of Dirac operators on 
the interval [0, 1]. Namely, we consider one-dimensional Dirac operators in the AKNS 
normal form [1] generated by the differential expressions 



with potentials 



CN ■ (1-1) Qi^) = 9i(a;)cri + q3{x)a3. 



Here 



(N 



^i=li ' i ' ^^-vo -1 



Q I are the Pauli matrices and qi and q2 are real-valued integrable functions that incorpo- 

rate various physical characteristics of the particle considered (its mass m, the scalar 
potential, electrostatic potential, and anomalous magnetic moment, see [Ml Ch. 4]). 
We assume that the potential is of bounded support and thus consider the Dirac op- 
erator on a bounded interval (which without loss of generality might be taken [0, 1]). 

The Dirac equation has widely been used in various areas of physics and mathematics 
starting from 1929, when P. Dirac suggested it to model the evolution of spin-i particles 
in the relativistic quantum mechanics f30]. In 1973 Ablowitz, Kaup, Newell, and 
Segur [1] observed that the Dirac equation is related to a nonlinear wave equation 
(the "modified Korteweg-de Vries equation", a member of the AKNS-ZS hierarchy, 
see [21 [33]) in the same manner as the Schrodinger equation is related to the KdV 
equations, and this discovery stimulated the increasing interest in direct and inverse 
problems for Dirac operators in both physical and mathematical literature. The topic 
was brought to life somewhat earlier; namely, in 1966, Gasymov and Levitan solved 
the inverse problems for Dirac operators on M4. by the spectral function [llj and by 
the scattering phase (TU]. Their investigations were continued and further developed 
in many directions. The reference list on the inverse spectral and scattering theory 
for Dirac operators has become very vast by now, and we only mention papers [SlEl 
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[iniES] that contain extended bibliography on the topic. The books by Levitan and 
Sargsjan [23] and by Thaller [3D] may serve as a good introduction into the (respectively 
mathematical and physical part of the) theory of Dirac operators. 

In what follows, we fix 9o,9i G [0,7r) and define the Dirac operator ^{9q,9i,Q) in 
the Hilbert space H := L2{0, 1) x ^2(0, 1) as the restriction of iq onto the domain 

dom^{9o,9i,Q) = {u = {ui,U2)' G H | ui,U2 G AC[0,1], 

ui{j) cos9j = U2{j) sm9j,j = 0, l}. 

It is well known that so defined operator ^{9o,9i,Q) is self-adjoint and has a simple 
discrete spectrum tending to ±00. In general, the spectrum of the operator ^{9q, 9i, Q) 
does not determine the potential Q. Indeed, there are only several analogues of the 
Ambartsumyan result for the Sturm-Liouville operators on a finite interval [1], where 
the unperturbed Dirac operator (i.e., with Q = 0) is the only one having the specific 
spectrum; see [TB ] [25 ] [3T| . To reconstruct the operator ^{9o, 9i, Q) in a generic situation 
some additional information is needed. For instance, the classical result of the inverse 
theory for the Dirac operators states that if we also know the spectrum of another op- 
erator ^(6*0, 9[, Q), with 9[ such that sin(6'^ — 9i) 7^ 0, then two such spectra determine 
uniquely the potential Q in the AKNS form (11. ip . 

For the sake of definiteness (and in order to avoid inessential technicalities), we 
choose the boundary conditions corresponding to = ^1 = and 9[ = it/ 2 and 
denote the operators ^{0,0, Q) and ^{0,n/2,Q) by ^i{Q) and ^2{Q) respectively. 
If Q G 1^2(0, 1) ® C^, then the eigenvalues (A„)„gz and {^n)n£Z of ^i{Q) and 2>2{Q) 
respectively can be enumerated so that they satisfy the interlacing condition 

(1.2) An < yUn < A„+i, n G Z, 

Tin + A^, 
7r(n + I) + 

where the remainders A^ and n'^ form ^2-sequences; see [5l[24]. Conversely, it was 
shown in [3] that any two sequence of real numbers that interlace and obey the asymp- 
totics (11. 3p are eigenvalues of the operators S^iiQ) and ^2(Q), for some real- valued 
potential Q in the AKNS form with L2-entries. 
We thus have a well-defined mapping 

(1.4) ((A„),(/^„)) 

from the spectral data ((An), (yUn)) of Dirac operators S'i{Q) and ^2{Q) to the corre- 
sponding potential Q in the AKNS form. A natural question arises, what properties 
this inverse spectral mapping has. 

Unfortunately, we are not aware of any work where this topic has been explicitly 
studied. However, there have been many papers discussing the analogous question for 
Sturm-Liouville operators with integrable potentials; see [6l[8l|27ll28] and the refer- 
ences cited in |T^ . (We observe in passing that the case of Sturm-Liouville operators 
in impedance form considered in [HIHlllHlET] is in a certain sense equivalent to the class 
of reduced Dirac operators with ^3 = 0; cf. [3].) Adapting as necessary the correspond- 
ing proofs for the Sturm-Liouville operators, one can show that the inverse spectral 
problem is locally stable, i.e., that the potential Q depends continuously on the spec- 
tral data ((A„), (fJ^n))- Here the metric on the set of the spectral data is induced by 



and the asymptotics 

A 

(1.3) 



n 
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that of £2 X ^2 on the elements ((A^), (/i^)) with the and /x^ defined by (II .Sp . More 
exactly (cf. [6],|27] for the Sturm-Liouville case), this local stability means that, for a 
fixed M > 0, there are positive e and L such that if ||(5i||h < ^5 ||Q2||e < and 
their spectral data ui := ((Ai,„), (/xi,™)) and 1^2 ■= ((A2,n), (/^2,n)) satisfy Hi^i -1/211 < e, 
then 

(1.5) ||Ql-g2||H<^||l^l-i^2||. 

Such a result cannot be considered satisfactory e.g. for the purpose of numerical 
reconstruction, as it refers to the norm of the potential Q to be recovered and thus 
specifies neither the allowed noise level e nor the Lipschitz constant L. Therefore it is 
desirable to have a uniform stability that asserts fll.Sp whenever the spectral data 1^1 
and U2 run through a bounded set ^ and with L only depending on 

Recently, such a uniform stability in the inverse spectral problem for Sturm-Liouville 
operators on [0,1] was established by Shkalikov and Savchuk [25] ■ They considered 
operators with real-valued potentials from the Sobolev spaces M^2^(0, 1) with s > —1; 
for negative s, such potentials are distributions. Their approach was based on a careful 
study of the nonlinear mapping of (11.40 and used extensively the implicit function 
theorem. In our work [19] we established analyticity and global stability of the inverse 
spectral mapping for s G [—1,0] using a different approach that generalizes a classical 
method due to Gelfand and Levitan [12] and Marchenko [26j and has been successfully 
apphed to reconstruction of Sturm-Liouville operators with singular potentials [2UII2T] . 
In [12], similar results were established for Sturm-Liouville operators in impedance 
form. 

The main aim of this paper (Theorem 12. ip is to prove analyticity and uniform con- 
tinuity of the inverse spectral mapping ((A„), (/U„)) Q for the class of the Dirac 
operators under consideration. In fact, we first establish these properties for the in- 
verse problem of reconstructing Q from the eigenvalues of (A„) and the corresponding 
norming constants (a^) introduced below (Theorem 12.20 . and then show that these 
norming constants depend analytically an Lipschitz continuously on two spectra, (A„) 
and (fin)- We use the approach to the inverse spectral problem for Dirac equation 
based on the Krein equation [5] and further develop the methods of [T^j. 

We mention that the methods of [3] could be used to treat the Dirac operators in 
AKNS form with real-valued potentials belonging to Lp(0, 1) ® with p G [l,oo). 
However, apart from some technicalities caused by more complicated properties of the 
Fourier transform in Lp{0, 1) for p ^ 2 the approach would remain the same and we 
decided to sacrifice the generality to simplicity of presentation. 

The paper is organised as follows. In the next section we introduce all the related 
objects, formulate the main results, and recall the approach to solution of the inverse 
spectral problem for Dirac operators in AKNS form based on the Gelfand-Levitan- 
Marchenko and Krein equations. In Section [31 we prove Theorem 12.21 and in Section H] 
study dependence of the spectrum (/i„) of the operator ^2{Q) on (A„) and (a„) and 
then derive Theorem l2.1[ Finally, three appendices discuss properties of special analytic 
mappings between Banach spaces and recall some facts on the special Banach algebra 
that are essentially used in the paper. 

Notations. Throughout the paper, we shall denote by L2(0, 1) the set of all 2 x 2 
potentials in the AKNS normal form; clearly, this space is unitarily equivalent to 
L2(0, 1) X 1^2(0, 1). Bold letters will usually denote vectors and functions with values 
in H, and ||/|| and ||f|| will stand for the L2(0, l)-norm of a scalar function / and 
H-norm of the vector-valued function f , respectively. 
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2. Preliminaries and main results 

In this section we state the main resuhs and recall the method of solution of the 
inverse spectral problem based on the Gelfand-Levitan-Marchenko [21] and Krein [3] 
equations. 

2.1. Spectral data. We denote by A„ and n e Z, the eigenvalues of the operators 
^i{Q) and ^2{Q) respectively and recall that these eigenvalues interlace, i.e., A„ < 
/^n < ^n+i for all n e Z, and satisfy the relations 

A„ = 7rn + P2n, /in = 7r(n + i) +P2n+1 

with some ^2-sequence (p„)- 

The equation £qu = Au, u = (mi,M2)*, subject to the initial conditions mi(0) = 0, 
^2(0) = 1 has the solution 



(2.1) s(x. A) = So(x, A) + / K{x,t)so{t, X) ds, 

Jo 

where s(x. A) := (si(x. A), 52(0;, A))*, So(x, A) := (sin Ax, cos Ax)* and K =: {kjiY- i^^ 
is the kernel of the transformation operator. Further, s(-,Ari) is an eigenfunction 
corresponding to the eigenvalue A„ of the operator &i{Q), and we call the number 
an := ||s(-, A„)||~^ the norming constant for this eigenvalue. It is known that 



where the sequence (a^)„gz belongs to £2 := ^2(^)- Moreover, the norming constants a„ 
can be determined from the spectra of the operators &i{Q) and Sl^iff) as follows. 

We set ^(A) := Si(l, A) and C(A) := ^2(1, A); due to (12. ip these are entire functions 
of exponential type 1 with zeros A„ and p„ respectively. The Hadamard canonical 
products of S and C are 



00 



, 00 

'An — Z r^/ \ ^.T TT t''n — Z 



(2.2) S{z) = {z- Ao)V.p. n , C{z) = V.p. J] 



nn ' ' ' ' J. J. jcffi _|_ l\ 

n=—oo ^ ^' 



(the prime denoting that the factor corresponding to n = is omitted), so that S and 
C are uniquely determined by their zeros. Then we have [3| 

(2.3) a, ^ 



where the dot denotes the derivative in z. 

2.2. Main results. We introduce the set J/' of data ((Ari)„gz, (A''n)nez) with the fol- 
lowing properties: 

• the sequences (A„) and (/i„) strictly interlace, i.e., A„ < /i„ < A„+i for all n G Z; 

• the sequence {pk)k&, with p2n := A„ - vrn and p2n+i := Pn - vr(n + i), belongs 
to £2. 

In this way every element u := ((A„), (/i„)) of ^ is identified with a sequence (p„) 
in £2 thus inducing a metric on 

According to [3] , every element of ^ gives the eigenvalue sequences of the operators 
^i(Q) and ^^{Q) corresponding to a unique real- valued AKNS potential Q G L2(0, 1) 
and, conversely, for every real- valued Q G L2(0, 1) of the form (11.11) the spectra of 
the corresponding Dirac operators form an element of ^ . When the AKNS poten- 
tial Q varies over a bounded subset of L2(0, 1), then the corresponding spectral data 
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((An), (/i„)) remain in a bounded subset of ,JV . Moreover, the Priifer angle technique 
yields then a positive h such that all the corresponding spectral data ((A„), (/U„)) are 
/i-separated, i.e., that /i^+i — A„ > /i and A„ — /x„ > /i for every n G N. Summarizing, 
we conclude that the uniform stability of the inverse spectral problem we would like to 
establish is only possible on the convex closed sets ^{h,r) of spectral data consisting 
of all elements of ^ that are /i-separated and satisfy ||(pn)||^2 ^ ^■ 
In these notations, the first main results of the paper reads as follows. 

Theorem 2.1. For every h G (0,7r/2) and r > 0, the inverse spectral mapping 

^{h,r) 3iy^Q e L2(0,l) 

is analytic and Lipschitz continuous. 



In fact, as in [19], we prove first the analyticity and Lipschitz continuity of the 
inverse spectral problem of reconstructing Q from the Dirichlet spectrum (A„) and the 
norming constants (a^) (see Theorem 12.21 below), and then derive Theorem 12.11 by 
showing that the norming constants depend analytically and Lipschitz continuously on 
the two spectra. 

More exactly, we denote by =Sf the family of strictly increasing sequences A := (A„) 
such that p2n '■= A„ — vrn forms an element of £2 and pull back the topology on =Sf 
from that of £2 by identifying such A with (p2n) ^ For h G (0, tt) and r > 0, we 
denote by Jtf^h, r) the closed convex subset of ^ consisting of sequences {Xn)n€Z with 
An+i — > h, and such that ||(p2n)i|£2 — ^- Next, we write ^ for the set of sequences 
a := (a„,)„gN of positive numbers (a„) such that the sequence (/3„) with /3„ := a„ — 1 
belongs to £2- This induces the topology of £2 on £/; we further consider closed subsets 
£/{h,r) of £^ consisting of all (a„) such that «„ > h for all n G N and ||(/9n)||^2 — 

It is known that, given an element (A, a) G =Sf x =2/, there is a unique Q G L2(0, 1) in 
the AKNS normal form such that A is the sequence of eigenvalues and o: the sequence 
of norming constants for the Dirac operator ^i{Q). Some further properties of the 
induced mapping are described in the following theorem. 

Theorem 2.2. For every h G (0, tt) and every positive h' , r, and r' , the inverse spectral 
mapping 

^(/i, r) X £^{h', r') 9 (A, a) ^ g G L2(0, 1) 
is analytic and Lipschitz continuous. 

2.3. Solution of the inverse spectral problem via the GLM and Krein equa- 
tions. The Gelfand-Levitan-Marchenko equation relates the spectral data for the op- 
erator &i{Q) (i.e., its eigenvalues and norming constants) with the transformation 
operator + To derive it, we start with the resolution of identity for ^i{Q), 

N 

= s-lim a„(-,s„)HS„, 

n=-N 

where s-lim stands for the limit in the strong operator topology and s„(x) := s(x, A„). 
Recalling that s„ = (J^ + ^)so,n with So,n(x) = So(x, A„), we get 



N 

2.4) J = + [s-lim V a„(-, So,„)eSo. J (^ + ^ 



n=l 
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The operator in the square brackets has the form J' + where ^ is an integral 
operator of Hilbert-Schmidt class with kernel 

(2.5) F{x,t):=\{E{^)-H{^)a,l 
and 

oo 

(2.6) H{s) ■= V.p. ^ (^a„e2^"'""^ - e^™^^ 

n=— oo 

Applying (j^ + ^*)~^ to both sides of (12.41) and rewriting the resulting relation in 
terms of the kernels K and F, we get the Gelfand-Levitan-Marchenko (GLM) equation 

(2.7) K{x,t) + F{x,t) + K{x,s)F{s,t)ds = 0, x>y. 

Jo 

If the potential Q is continuous, then the algorithm for reconstructing Q from the 
spectral data ((A„), (/U„)) proceeds as follows, cf. [23]. We first calculate numbers a„ 
via (12. 3p . then construct a matrix-function H of (12. 6p . form the kernel F of (12. 5p . solve 
the GLM equation (12. 7p . and finally set 

(2.8) Q{x) = K{x,x)-(72 a2K{x,x). 

i i 

If the potential Q is only assumed to be in L2(0, 1), then the kernel K has the 
property that K{x,-) and K{-,x) are square-integrable matrix-functions depending 
continuously in ^2(0, 1) ® C'* on a; G [0, 1]. In particular, K may not have well defined 
restriction K{x, x) on the diagonal and thus the above-described method cannot be 
used. 

To overcome this obstacle, we consider, along with the GLM equation, the Krein 
equation 



(2.9) R{x,t) + H{x -t) + R{x,s)H{s-t)ds = 0, < t < x < 1. 

Jo 

It is easily seen that if i? is a solution to this equation, then 

K{x, t) := i [R{x, ^) - R{x, ^)as] 

solves the GLM equation, and in this sense the Krein equation is more general. The 
kernel R generates another transformation operator in the sense that 

r 

(2.10) s(x. A) = So(x, A) + / R{x,x — t)sQ{x — 2t, \) ds. 

Jo 

Moreover, as soon as the related convolution operator 

(2.11) J^u{x) := [ H{x-t)u{t)dt 

Jo 

satisfies the inequality + Jif > in ^2(0, 1) x C^, the Krein equation can be shown 
to possess a unique solution R and this solution has the property that R{-,t) belongs 
to 1/2(0, 1) ® C"^ and depends continuously therein on t G [0, 1]. Now, rewriting the 
relation (12. 8p in terms of R, we find that 

(2.12) Q{x) = -R{x,0)ai. 

This formula defines Q as an L2(0, 1) ® C^- valued function and will be the basis of our 
reconstruction algorithm. 
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3. Stability of the inverse spectral problem: norming constants 

In this section, we prove Theorem 12.21 on uniform stabihty of the inverse spectral 
problem of reconstructing the potential Q of a Dirac operator ^i{Q) from its spectrum 
and norming constants. 

We shall study the correspondence between the data (A, a) G J^{h, r) x [h', r') 
and the potentials Q in the AKNS form of the Dirac operator ^i{Q) through the chain 

{X, a.) H R Q , 

in which H is the matrix- valued function of fl2.6p . R is the matrix-function solving the 
Krein equation (12. 9p . and, finally, Q is given by fl2.12p . 

Lemma 3.1. The mapping 

^{h, r) X £^{h', r') 3{X,a)^H e L2{-1, 1) ® 

is analytic and Lipschitz continuous. 

Proof. Since the matrix B has simple eigenvalues ±i, it is unitarily equivalent to 
the matrix diag{i, — i}, and thus H is unitarily equivalent to the matrix-function 
diag{/i(s), s)}, where 

oo 

(3.1) h{s) := V.p. J2 («ne^^"'' " e^™')- 

n=— oo 

We thus need to prove that the mappings 

^{h, r) X £/{h', r') 3 (A, a) ^ h{±s) G L2{-1, 1) 

are analytic and Lipschitz continuous. It clearly suffices to consider only the case 
of h{s)] also, since the a„ and A„ are real, we get h{—s) = h{s), and only the restriction 
of h onto (0, 1) need to be studied. 
We have h = hx + hx^a, where 

oo oo 

hx{s) := V.p. J2 [e'"'"" - h^As) ■■= V.p. f^ne^'"'". 

n=— oo n=— oo 

We recall that the numbers /3„ := a„ — 1 and p2n = A„ — nn form sequences in £2 that 
induce the topology of ^ and The series 



V.p. P2ne'™" and V.p. ^ /^^e^™^ 

neZ 



converge in ^2(0, 1), and we denote by f\ and ga respectively the corresponding sums. 
Since the sequence of functions forms an orthonormal basis of -^2(0, 1), the 

mappings (p2n) ^ fx and {/3n) ga are unitary isomorphisms between £2(^) and 
L2(0,l) given by the inverse Fourier transform. Therefore we have hx = g{fx) and 
hx,a = ^{f\,ga) with the mappings g{f) and ^{f,g) introduced in Appendix [Bl The 
claim now follows from the properties of g{f) and ^{f,g) established in Lemmata lB.il 
and IB. 21 respectively. □ 

Solubility of the Krein equation crucially relies on the following property of the 
convolution operator Jif. 
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Lemma 3.2. For every h G (0, vr), /i' G (0, 1), and positive r and r' , there exists e > 
with the following property: if (A, ck) is an arbitrary element of J^{h,r) x s^{hi,ri) 
and H is the function of (12.61) . then for the corresponding convolution operator we 
have + J^> eJ . 

Proof. Denote by ^± the convolution operators in L2(0, 1) constructed as in (12. lip 
but for the scalar functions h{±s) of (13. ip . Transforming H to the diagonal form as 
in the proof of the previous lemma, we see that it suffices to prove that there is £ > 
independent of (A, a) such that / + J^± > el. We shall only treat the case of the 
operator I + J^+, as the other case is completely analogous. 

Since the system (e^™'*)^^^ is an orthonormal basis of -^2(0, 1), we find that 



k 

((/ + jn)/,/) = (/,/)+ lim Yl [«n|(/,e^^-^)p-|(/,e2— )p] 

A;->-oo 

n=—k 

00 

= V.p. J2 «ni(/,e^^-^)r. 



The inclusion a G s^{h',r') implies that > h' for all n G Z. It follows from the 
results of [131 Ch. VI], |32l Ch. 4] that the system := (e^^"'')„g^ is a Riesz basis 
of Iv2(0, 1); moreover, there exists m = m{h,r) > that gives a lower bound of S'x for 
every A G ^{h,r) [T7]. Therefore, 

00 CO 

((J + ^+)/,/)=V.p. c^n\{f,e''"n\'>h'V.p. Y \U\e''n\'>h'm\\ft. 

n=—oo n=— 00 

The proof is complete. □ 

To study solvability of the Krein equation (12. 9p . we shall regard it as a relation 
between the corresponding integral operators. To this end we recall several notions 
that will be used. By definition, the ideal ©2 = ©2 (EI) of Hilbert-Schmidt operators 
consists of integral operators ^ in H whose kernels T satisfy the inequality 



7' 

Jo 



ti (T{x,y)T*{y,x)) dxdy < 00; 



here tr stands for the trace of a matrix or a trace-class operator. Notice that if T 
(^*i)L=i. then 



tr TT* = Y 1^' 



2 

2 



i.e., the operators in ©2 have kernels whose entries are square integrable in 17 = (0, 1) x 
(0, 1). The linear set ©2 becomes a Hilbert space under the scalar product 

{^,W)2:=tT{^W*)= [ [ ti {X{x,y)Y*{y,x)) dxdy 

Jo Jo 



UNIFORM STABILITY 



As noticed above, the matrix- valued function H satisfies the relation H*{s) = H{—s) 
and H{s) is unitarily equivalent to the diagonal matrix diag{h{s), h{—s)}. Therefore, 

ti H{x — y)H*{x — y) dx = / \h{x — y)]"^ dx + / \h{y — x)]"^ dx 



< 2 / \h{s)\'ds 



-1 



T+{x,y) 



so that the convolution operator Jif belongs to ©2 and its norm satisfies ||^||| = 

Denote by the subspace of 62 consisting of all Hilbert-Schmidt operators with 
lower-triangular kernels. In other words, ^ G ©2 belongs to ©^ if the kernel T of 3^ 
satisfies T{x,y) = for < x < y < 1. For an arbitrary G ©2 with kernel T the 
cut-off T"*" of T given by 

T(s, y) for X > y 
for X < y 

generates an operator G ©^, and the corresponding mapping : 3^ H- 
turns out to be an orthoprojector in ©2 onto ©^, i.e. (^'^)^ = and {^~^ 3^, '3^)2 = 
{^,^+W)2 for any e ©2; see details in Ch. I.IO]. 

With these notations, the Krein equation (12.91) for the kernels R and H can be recast 
as the relation between the corresponding Hilbert-Schmidt integral operators 

(3.2) ^ + + i^je) = 

or 

where is the linear operator in ©2 defined by = ^'^{'3^ ^) and 3 is the 

identity operator in ©2. Therefore solvability of the Krein equation and continuity of 
its solutions on is strongly connected with the properties of the operator 

Lemma 3.3. For every ^ G ^{T-L), the operator is bounded in ©2. Moreover, 
for every Jif from the set 

Sj:={jif = ^(A, a) I (A, a) G ^{h, r) x £^{h', r')} C ©2 

the operator 3 + is invertible in ^(©^) and the inverse {3 + depends ana- 
lytically and Lipschitz continuously on M' & in the topology of &2- 

Proof. Boundedness of is a straightforward consequence of the inequality 

cf. [HI Ch. 3]. Assume next that + ^ > eJ^ inM; then for ^ G ©^ we find that 

{{3 + ^J;)^, = + (^^, = tr (^ + JT)^*) . 

Since ^ ( -|- > and the trace is a monotone functional, we get 

i.e., 3^^\-> e3 in ©+. 

Applying now Lemma 13. 2[ we conclude that for every Jif G it holds 3 + > e3 
with e of that lemma depending only on h, h', r, and r'; therefore, 3 + ^~^ is boundedly 
invertible in ^(©^) and 
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Since depends linearly on J^, it follows that the mapping M' i— )■ (J + 

from ©2 into is analytic and Lipschitz continuous on the set i^. The proof is 

complete. □ 

Corollary 3.4. For every M' G io, the Krein equation ( \2.9\\ has a unique solution 

M := -{3 + G et; 

moreover, ^ depends analytically and Lipschitz continuously in on E Sj G &2- 

It follows that the kernel R{x,t) of M is square integrable in the domain Q and 
depends analytically and Lipschitz continuously in L2{fl) ® on J^. However, we 
need to know that R{x, 0) is well defined and belongs to L2{0, 1) Cg) C^. 

To this end we use the Krein equation to show that 

R(x, t) = -H{x -t)- f R{x, s)H{s - t) ds 

Jo 

as a function of x depends continuously in L2{VL) x on t G [0, 1]. Since the shift 
/(■) ^ f{'~^) is a continuous operation in L2(M), H{-—t) enjoys the required property. 
Denote by Rij and Hij, i,j = 1, 2, the entries of the matrix-valued functions R and H 
respectively. They belong to L2{fl) since the corresponding integral operators ^ and 
Jif are in 62; therefore, 



Rij{x, s)Hjk{s — t) ds 







2 



dx 



(3.3) < 



dx \Rij{x, s)]"^ ds / \Hjk{s — t)]"^ ds 
Jo Jo Jo 



•1 

< / \Hjk{s)\'^ ds j j \Rij{x, s)\ ds dx < 00. 



Thus the function 

(3.4) [ Rij{x,s)Hjk{s-t)ds 

Jo 

of the variable x G [0, 1] belongs to ^2(0, 1); moreover, continuity of the shifts Hjk{- —t) 
and estimate ( 13. 3p show that function f l3.4p depends continuously in L2(0,l) on t G 
[0, 1]. Since the entries of the matrix R{x, s)H{s—t) ds are sums of functions ( 13. 4p for 
corresponding indices i, j and k, we conclude that R{-,t) indeed depends continuously 
in L2(0, 1) C"^ on t G [0,1]. In particular, i?( ■ , 0) is a well-defined function in 
L2(0,1)®C^ 

Proof of Theorem \2.S[ The above arguments show that it is legitimate to take t = in 
the Krein equation; thus we get that 



R(x, 0) = -H{x) - I R{x, s)H{s) ds. 
Jo 



The right-hand side is a bilinear expression in H and R. In view of the analytic depen- 
dence of 3? onto Jif stated in Corollary 13.41 and estimates fl3.3p . this yields analyticity 
and Lipschitz continuity of R{x, 0) on if G L2(— 1, 1) (S) C"^. By (12.121) the function Q 
enjoys the same properties, and the proof is complete. □ 
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4. Reconstruction from two spectra 

We recall that the norming constants a„ for the Dirac operator ^i{Q) can be deter- 
mined from the spectra (A„) and of ^i{Q) and ^2{Q) by the formula 

1 

5(A„)C(A„)' 

where the entire functions S and C are given by the canonical products (12 ■2p over A„ 
and iJ,n respectively. This induces a mapping i/ i— )■ o: from the spectral data u : = 
((A„), (/x„)) e into the norming constants a := (an) € In this section, we shall 
establish Theorem 12. II by proving the following result. 

Theorem 4.1. For every h G (0,7r/2) and r > 0, the mapping 

(4.1) J^{h, r) 3 v a. e 

is analytic and Lipschitz continuous; moreover, there exist positive constants h' and r' 
such that the range of this mapping belongs to £^{h',r'). 

By definition, ^ consists of elements of the commutative unital Banach algebra A 
introduced in Appendix O We observe that the metrics on agrees with the norm 
of A, and thus the results of Appendix [C] yield the following statement. 

Proposition 4.2. For every positive h and r, the set ^{h,r) consists of invertible 
elements of A. Moreover, the mapping a i— ck^^ is analytic and Lipschitz continuous 
in A on ^ {h, r), and its range lies in .g/ ((1 + r)^^ ,rh^^^ . 

In view of Proposition 14. 2^ it suffices to prove Theorem 14. II with ex replaced by (X~^. 
The elements of the sequence ck"^ are = S'(A„)C(A„). We shall show that the 
sequences 

7 := ((-1)"^(A„))„^,, 6 := ((-1)"C(A„,))„^, 
form elements of Thus Theorem 14.11 will be proved if we show that the mappings 

(4.2) ^{h, r) 9 J/ i-> 7 G ^{h, r) 3 v ^ 5 e 

enjoy the properties required for the mapping (14. ip . 

To begin with, integral representation (I2.10p of the solution s(-. A) yields the formulae 

(4.3) S{X) = smX+ [ r^{s)e'^'^^-^'Us, 

Jo 

(4.4) C(A) = cosA+ / r2{s)e'^^''^'Us 

Jo 

for the functions 5* and C. Here vi and r2 are some (uniquely defined) functions 
in -^2(0, 1) formed from the entries of the matrix-valued function -R(l, 1 — x). Therefore 
both S and C can be recast in the form 

cosA+ / g{s)e'^^^-^'Us 
Jo 

with g{s) = i(l — 2s)ri(s) for 5* and g{s) = r2(s) for C. The sequences 7 and <5 have 
therefore similar structures; namely, their n-th element equals 

(4.5) cosp2n + (-1)" [\{s)e'^"^^-^'^ ds 

Jo 

for respective g; here, as usual, p2n := A„ — vrn. 
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Clearly, the mapping (p2n) ^ (cosp2n — 1) is analytic in £2- Its Lipschitz continuity 
follows from the inequality | cosa; — cos?/| < |x — ?/|; also, the inequality 1 — cosx < 
yields the estimate ||(cosp2n - 1)|| < ||(P2n)|P/2. 

Let fx denote the function in 1/2(0, 1) such that f\{n) = p2n and set 

(4.6) K:={-ir [\{s)e''-'^'-''Us. 

Jo 

Then the function h given by the Fourier series Ylnez hn^'^™''^ coincides with the func- 
tion \E'(/a,(7) of Lemma lB.31 It follows from that lemma that the sequence (/i„)nez 
depends analytically and boundedly Lipschitz continuously in £2 on fx and g. We 
prove in the lemma below that the function g depends in the same manner on 1/ = 
(A,M)e^(/i,r). 

Lemma 4.3. The mappings 

^{Kr) 3 (A,/x) ^ r, e L2(0,l), j = 1,2, 
are analytic and Lipschitz continuous. 

Proof. Since both mappings can be treated similarly, we only discuss the first one. By 
definition, we have S{Xn) = 0, and thus the numbers A„ = nn + p2n, n E Z, are zeros of 
the entire function S of 04.31) . The required properties of the mapping u ^ ri follow 
now from the results of Appendix [XI □ 

The above reasoning justifies the inclusion ol^^ G ^ as well as analyticity and 
Lipschitz continuity of the mappings of (14.21) . It remains to prove that there exist 
positive h' and r' such that, for every u G ^{h,r), the corresponding elements 7 and 
5 belong to £^{h\r'). 

Existence of such an r' follows from the uniform estimates of the £2-norms of the 
sequences (cosp2„ — 1) and (/i„); see (14. 5 p and (14.61) . Indeed, in view of Lemma fB.SI 
the function h = Xlnez remains in the bounded subset of L2(0,l) when fx 

and g vary over bounded subsets of -^2(0, 1), and the latter is the case when u runs 
over ^{h,r) by the definition of the functions fx and g and Lemma [4. 3[ 

Next, in view of formula (12.21) and the interlacing property of A„ and p^, the numbers 
7„ = (— 1)"S'(A„) and 5„ = (— 1)"C(A„) are all of the same sign and thus are all positive 
in view of the asymptotic relation (14.51) . The uniform positivity of 7„ and 5„ (and thus 
existence of a positive h' such that = 7„5„ > h'), immediately follows from the 

lemma below. A similar statement is established in jT3|; we give a complete proof 
below for the sake of completeness. 

Lemma 4.4. For every h G (0, 7r/2) and r > we have 

sup sup log [^(An)! < 00, sup sup log |C(A„)| < 00, 

(A,/x) neN (A,/i) nGN 

where S and C are constructed via (12.21) from sequences X and fi, and the suprema are 
taken over (A, /x) G cA^{h, r). 

Proof. We assume first that n ^ 0. By (12. 2p . we hav^ 

A/ \ N _ Ao — A„ T-r Afc — A„ 
0[An) — I I ; • 

TTU vr/e 



In what follows, all summations and multiplications over the index set Z will be taken in the 
principal value sense and the symbol V.p. will be omitted. 
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Dividing both sides by 



dsinz 



cos 7m 



dz 



n 



7rk — Tin 



we conclude that 



i^(A„)i= n 



keZ, k^O,n 



Afc — A„ 
7T{k — n) 



Ilk 



and observe that this formula also holds for n = 0. Set (recall that p2k '■= — T^k) 

^k ~ ^ _ P2k — P2n 

0'k,n '■— — ^ J- 



TT(k — n) 



TT 



(A; — n) 



ii k ^ n and a„„ := 0; then [^(An,)! = Ilfcezl-'- + '^k,n)- Since the sequence (A„) is 
2/i-separated for every (A, /x) G r), we have 1 + au^n > 2/i/7r for all integer k 

and n. Therefore, with 



K := max 

x>-l+2h/'K 



log(l + x) — X 



x^ 



< oo, 



we get the estimate 
(4.7) 



log JJ(1 + afc,n) <\^a'k,n +K^: 



'k.ni 



k& 



feez 



kel 



provided the two series converge. 
Clearly, 



k^n 



k — n 



and thus 



fcez 



P2k 



TT k — n 

k^n 



< 



by the Cauchy-Bunyakovski-Schwarz inequality (recall that XlfcgzPifc — '"^ by the 
definition of the set ^{h, r) and Xlfc^nl^ ~ '"■)~^ ~ Next, the inequality 



2 < ^Plk 



+ 



- (A; - n)2 ^ (fc - n) 



for k ^ n yields 



2^2 



k& 



k^n 



(n - A;)^ 



It follows from f l47fl) that 

llog JJ(1 + afc,0 < ( VSr + 4i^7rV2)/3, 



fcGZ 



where the constant K only depends on h. 
Similarly, we find that 

Pk — \. 



ic(A..)iHn£^ =n 



pk ~ An 



and then mimic the above reasoning to establish the other uniform bound. The lemma 
is proved. □ 
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Combination of the above results constitutes a complete proof of Theorem 14.11 and 
together with Theorem 12.21 establishes Theorem 12. 1[ 

Acknowledgements. The author thanks A. A. Shkalikov and Ya. V. Mykytyuk for 
stimulating discussions. The research was partially supported by the Alexander von 
Humboldt Foundation and was partially carried out during the visit to the Institute for 
Applied Mathematics of Bonn University, whose hospitality is warmly acknowledged. 

Appendix A. Analyticity of some related mappings 



Here we give a brief account on the results of [22] as well as some extensions that 
are used to prove Lemma [4.3[ It was shown in [22j that for every / G i^2(0, 1) there 
exists a unique function g G i^2(0, 1) such that all zeros (counting multiplicities) of the 
entire function 

(A.l) Gg{z):=siYiz+ ! gity'^^-^'Ut 

Jo 

are given by the numbers im + f{n), n G Z. Such pairs of / and g in fact satisfy the 
relation 

(A.2) H{f, g) := s{f) +g + J2 fc, = 0^ 

k=l 

here 

^^^^ '"SZ (2k + 1)1 ' 

k=o ^ ' 

f^^^ is the fc-fold convolution of / with itself, and M is the operator of multiplication 
by i(l — 2x). The function H is analytic from ^2(0, 1) x ^2(0, 1) into ^2(0, 1), and its 
partial derivatives dfH{f, g) and dgH{f, g) are given by 

(A.3) dfHif, g){h,) = (c(/) + J2 ^ ) * hi, 



k=l 

00 



{k-l)\ 



(A.4) dgH{f,g){h,) = h, + J2 ^, 

k=l 

with 



=(/) - E 



k £<2k> 



Using the implicit function theorem, it was shown that the induced mapping f g is 
analytic. In order to establish its Lipschitz continuity, we shall study the above partial 
derivatives in more detail. 

Namely, we assume that / G ^2(0, 1) is such that the corresponding sequence A = 
(An)ngN with An := vrn + /(n) belongs to ^{h,r). Set Sx to be the canonical product 
of fl2.2p : then Sx can also be represented as (lA.ip . Direct calculations show that the 
n-th Fourier coefficient of the function of (IA.3P is equal to 

(-l)"/ii(n)[cosA„+ [ i{l-2t)g{t)e'^-^^-^'Ut\ hi{n)Sx{Xn). 



By Lemma 14.41 there are positive numbers Ki and K2 such that 

Kl < \Sx{Xn)\ < K2 
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for all A G =Sf (/i, r) and all n & Z. Therefore the partial derivative dfH{f, g) is 
a bounded and boundedly invertible operator in ^2(0, 1); moreover, for every fixed 
h> and r > 0, the norms of dfH{f,g) and their inverses are uniformly bounded for 
/ G i^2(0, 1) generating the sequences A in the set ^{h,r). 

Similarly, the n-th Fourier coefficient of the function of ( 1A.4P is equal to 



"1 

-1)" / h2it)e'^-^^-^'Ut. 







By the results of [T7], there exist positive M and m such that, for all A G ^{h,r), 
the sequences (e''^"'-^~^^'')„gz form Riesz bases of -^2(0, 1) of upper bound M and lower 
bound m. Therefore the operator Hg := dgH{f,g), 



2nnix 



is bounded and boundedly invertible in L2(0,l), with \\Hg\\ < M^/^ and H-f^^ ^|| < 

We now use the implicit mapping theorem to conclude that the mapping f ^ g 
is analytic in 1/2(0, 1). The uniform bounds on the inverses of the partial derivatives 
dfH{f,g) and dgH{f,g) established above imply that, for every h> and r > 0, this 
mapping is Lipschitz continuous on the set of functions / G L2(0, 1) generating the 
sequences A G ^{h, r). 

Appendix B. Some auxiliary results 

We recall that the convolution f * g oi two functions in L2(0, 1) is a function in 
L2(0, 1) given by 

{f*9){x) := I f{x-t)g{t)dt, 
Jo 

where / is extended to (—1,0) as a periodic function with period 1. The (discrete) 
Fourier transform / of / G 1^2 (0, 1) is a function over Z given by 

f{n) := / f{t)e~^^'''Ut. 







It is well known that the Fourier transform is a unitary mapping from ^2(0, 1) to ^'2(^1) 
and that / * g{n) = f{n)g{n)] as a result, we have the inequality 

ll/*^7ll< 11/1111^711 

for all /,(7GL2(0,1). 

Lemma B.l. For a function f G ^^2(0, 1), set 



2nnix 



g{f){x):=Y.p. J2 [e'^"^"^" - l]e 

n=— 00 

Then the series determines a function in -^2(0, 1), and the mapping 
is analytic and locally Lipschitz continuous on bounded subsets. 
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Proof. We start with observing that the series Ylnez f '^i^)^'^^™'^ Fourier series 

for the function f^''\ the /c-fold convolution of / with itself, and that ||/^'^^|| < " ^"'^ 
Developing e-^*-"-'''* into the Taylor series, we find that 



g{f) = v.p. y: [e 



oo 



n=—oo k=l 
\k 



fc=l 



oo 



f{k) 
k\ ■' ■ 

k=l 

The change of the summation order in the second equality above is justified by the fact 
that, for k > 1, the summands in the double series are dominated by C^p{n)/k\ with 
C := 2 max„gz{|/(r2)|} + 1. Therefore the double series over the index set {{n, k) \ n E 
Z,k > 1} converges absolutely and the Fubini theorem applies. This formula represents 
g{f) as an absolutely convergent series (which is a Taylor series expansion of g{f) in 
the variable /) and thus proves the analyticity in ^2(0, 1) of the mapping / 1— g{f). 
Lipschitz continuity of that mapping on bounded sets follows from the estimate 



\9ih)-9im = \\j:'-^uf-f?'] 

k=l 

00 



fc=l ^ ^' 

which is valid as soon as the L2-norms of /i and /2 are not greater than r. The proof 
is complete. □ 

Lemma B.2. For f and g in -^2(0, 1), set 

00 

*(/,^):=V.p. E ^Hexp{2[7rn + /(n)]is}. 

n=— 00 

Then the function $(/, (7) belongs to L2{0, 1) and the mapping 

$ : 1.2(0, 1) X 1.2(0,1)^^2(0,1) 
is analytic and Lipschitz continuous on bounded subsets. 

Proof. Transformations similar to those used in the proof of the above lemma show 
that 

^(/'^) = E^[/^'^^*^]- 

k=l 

The mapping $ is linear (and thus analytic) in g, and its analyticity in / as well as 
Lipschitz continuity on bounded subsets is established in the same manner as for the 
mapping g{f) of Lemma [B.ll □ 

Lemma B.3. For f and g in 1.2(0, 1), set 

^{f,g) := V.p.5^(-1)" /^(t)exp{[7rn + /(n)]i(l-2t)}cite2™- 
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Then the function '^{f^g) belongs to -^2(0, 1) and the mapping 

^ : ^2(0,1) X L2(0,l) ^ ^2(0,1) 
is analytic and Lipschitz continuous on hounded subsets. 
Proof. The coefficient of e^'^'"^ in the above series for \I/ can be written as 

(B.l) /^(t)exp{i(l - 2t)/H}e-2™*rft 

Jo 

and gives the n-th Fourier coefficient of the function h := Yl'k'=o^k/k\, with Hq := g, 
hk := /^'^^ * M'^g for k > 1, and M being the operator of multiphcation by the 
function i(l — 2t). In other words, we have (7) = h. Since ||/ * g\\ < \\f\\\\g\\ for 
every / and g in ^2(0, 1), the functions belong to -^2(0, 1) and their norms there 
obey the estimate 

\\hk\\ < ll/f l|Mf ll^ll, 
with ||M|| denoting the norm of the operator M. Thus the series for h converges 
absolutely and, since every hk is a multi-linear function of / and g, the mapping \E' 
is analytic. Its Lipschitz continuity on bounded subsets is established in the usual 
manner, and the proof is complete. □ 

Appendix C. Banach algebras 

The space £2 = £i^) is a commutative Banach algebra under the pointwise multi- 
plication [xn) ■ iVn) = " Vn) ■ Its uuital extcusion A consists of elements of £00 of 
the form al + x with a G C, the unity 1 G £00 having all its elements equal to 1, and 
X = (xn) G £2- The norm in A is given by 

||al + x||a = \a\ + ||x||, 

and al + X is invertible in A if and only if a 7^ and a + x„ 7^ for all n G Z; in this 
case the inverse is equal to a~^l + y, where y = with y„ := —Xn/a{a + Since 
under the above assumptions we have inf„ |a + x„| > 0, we see that y indeed belongs 
to £2] moreover, 

||(al + x)"^|U < krVl + ||x||/inf |a + x„|). 

n 

The mapping x 1— x~^ is analytic on the open set of all invertible elements of A; in 
addition, it is Lipschitz continuous on the sets 

=5^e := {al -|- X I |a| > 5, inf \a + x„| > e}. 

n 
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